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ABSTRACT 

We report on the discovery of GJ 3470 b, a transiting hot Uranus of mass m p = 14.0 +1.8 M e , radius R p = 4.2 + 0.6 R© and 
period P = 3.3371 + 0.0002 day. Its host star is a nearby (d=25.2+2.9pc) M1.5 dwarf of mass M* = 0.54 + 0.07 M Q and radius 
R+ = 0.50 + 0.06 R Q . The detection originates from a radial-velocity campaign with Harps that focused on the search for short-period 
planets orbiting M dwarfs. Once the planet was discovered and the transit- search window narrowed to about 10% of an orbital period, 
a photometric search started with Trappist and quickly detected the ingress of the planet. Additional observations with Trappist, 
EulerCam and Nites definitely confirmed the transiting nature of GJ3470b and allow for the determination of its true mass and 
radius.The star's visible or infrared brightness (V mag = 12.3, K™^ = 8.0), together with a large eclipse depth D = 0.57 + 0.05%, ranks 
GJ 3470 b among the most favorable planets for follow-up characterizations. 

Key words, stars: individual: GJ 3470 - stars: planetary systems - stars: late-type - technique: radial-velocity, transit 



1. Introduction 

Planets in transit are being detected by the thousands thanks 
to the Kepler mission (Borucki et al. 2011), with some be- 
ing smaller than Earth (Fressin et al. 2012; Muirhead et al. 
2012) and some lying in their star's habitable zone (Batalha 
et al. 2012). Kepler's detections remain however around too faint 
stars, and/or with too shallow transit depths, to perform transmis- 
sion and occultation spectroscopy of their atmosphere. Actually, 
most Kepler host stars are too faint to ascertain the planetary 
nature of most detections or to measure their mass with current 
spectrographs. And therefore, observations that aim to charac- 
terize the atmosphere of exoplanets kept focused on the planets 
that transit bright, nearby and small stars which, so far, have been 
detected with ground-based instruments. 

In particular, searches driven by radial- velocity (RV) obser- 
vations educate the photometric observations by providing tar- 
gets known to harbor planets, together with a time window to 
search for their possible transits. Targeting the small M dwarfs, 
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the photometric follow-up of every planet detection is easily 
done from the ground with a small aperture telescope (e.g. Gillon 
et al. 2007b; Nutzman & Charbonneau 2008). For instance, the 
transit of a Jupiter- (resp. Neptune-) like planet produces a 4-% 
(resp. half-%) drop in flux when crossing a 0.5-R© star. And, 
since the occurrence rate of Jupiter and Neptune-like planets 
with short orbital periods is low (below a few percent - Bonfils 
et al. 2011b), their a priori RV detection and ephemeris enhance 
the transit-discovery power of a small telescope by a factor of a 
few hundred. Although smaller size planets are more frequent, 
their ground based detection is more difficult. One must then fo- 
cus on the smallest M dwarfs (e.g. Charbonneau et al. 2009) or 
make use of a space-borne observatory (e.g. Bonfils et al. 201 lb; 
Demory etal. 2011). 



This paper reports on the detection of GJ 3470 b a transit- 
ing hot Uranus detected in the framework of our radial- velocity 
search with Harps and systematic photometric follow-up. Its 
outline is as follows : we present the stellar properties of GJ 3470 
in Sect. 2, the radial- velocity detection in Sect. 3 and the photo- 
metric detection in Sect. 4. Next, we derive the posterior proba- 
bility of stellar and planetary parameters for the system in Sect 5 
before presenting our conclusions in Sect. 6. 
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Table 1. Observed and inferred stellar parameters for GJ 3470. 
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Fig. 1. Upper panel: RV time series of GJ3470. Bottom panel: 
Periodogram of GJ 3470 RVs. The horizontal lines shows different level 
of false-alarm probabilities. 



2. The properties of GJ 3470 

GJ 3470 (aka LP 424-4, 2MASS J07590587+ 1523294, 
NLTT 18739) is a M1.5 dwarf (Reid et al. 1997) seen in the 
Cancri constellation (a = 07 / *59 m 06*, 6 = + 15°23'30") at a 
distance d = 25.2 ± 2.9 pc known via photometric calibration 
(Lepine & Shara 2005). 

The Simbad astronomical database also identifies GJ3470 
as Melotte 25 EGG 29, which classifies the star as a member 
of the Hyades cluster (Melotte 25) and could constrain its age. 
The denomination originates from Eggen (1990) who assigned 
GJ 3470 to the Hyades supercluster on the basis of its proper 
motion. However, Eggen also found that, to be a member of 
the Hyades supercluster, GJ 3470 would have a radial velocity 
of about +41.1 km/s. For few other stars in his sample, Eggen 
could compare his radial velocity predictions with literature val- 
ues and found good matches, most often with differences not 
exceeding few hundred m/s. For G J 3470 however, no literature 
radial velocity was available at the time and, today, we measure a 
very different value of +26.5 km/s. Thus, we could not consider 
GJ3470 as a member of the Hyades supercluster and sought 
clues from other age proxies. On the one hand, we looked at 
GJ3470's Ha (6562. 808A) and vsinz. The spectral line is seen 
in absorption and the projected rotational velocity is low (< 2 
km/s), both indicative of a mature star older than - 300 Myr. On 
the other hand, we transformed GJ3470's proper motion (Lepine 
2005) and systemic radial velocity (this paper) into galactic ve- 
locities (U=+30, V=-12, W=-l() km/s). We found they match 
the kinematic population of the young disk (Leggett 1992), what 
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0.38, 0.40 





Note that the stellar mass and radii are also constrained by the data 
of this paper and that posterior values are given in Table 2. R97: Reid 
et al. (1997); Lep05: Lepine & Shara (2005); R01: Rapaport et al. 
(2001); Cu03: Cutri et al. (2003); LegOl: Leggett et al. (2001); D00: 
Delfosse et al. (2000); B98: Baraffe et al. (1998); H97: Haywood et al. 
(1997); C100: Claret (2000) 



suggests an age of < 3 Gyr (Haywood et al. 1997) as well as an 
approximately solar metallicity. 

Rapaport et al. (2001) give an apparent brightness V mag = 
12.27 + 0.02 in an approximately Tycho band, which is com- 
patible with the V mag ~ 12.44 estimate from Lepine (2005).We 
used the photometric distance to convert its infrared photometry 
K = 7.989 + 0.023 (Cutri et al. 2003) into an absolute magnitude 
M™ ag = 5.98 + 0.58 and, after bolometric correction (BC™ g = 
2.59 - Leggett et al. 2001), into a L* = 0.029 + 0.002 L lumi- 
nosity. The K-band mass-luminosity relation of Delfosse et al. 
(2000) gives a 0.45 + 0.11 M© mass. Also, we use the theoreti- 
cal mass-radius relation of Baraffe et al. (1998), in good agree- 
ment with interferometric measurements (Demory et al. 2009). 
Assuming an age greater than 300 Myr and solar metallicity, we 
converted M+ into 7?* = 0.42 + 0.10. The same models also pro- 
vides mass-temperature relation and hence r ef f = 3600 + 200 K. 
To account for limb-darkening when modeling eclipses, we used 
the quadratic coefficients of Claret (2000) for T eff = 3500 K, 
logg = 4.5 and [Fe/H]=0 in the z-band (u a = 0.40, u b = 0.19) 
and V-band (u a = 0.38, u b = 0.40) filters. 



3. Spectroscopic detection 

We observed GJ 3470 with the Harps spectrograph, the state-of- 
the-art velocimeter fiber- fed by the ESO/3.6-m telescope (Mayor 
et al. 2003; Pepe et al. 2004). Our settings remained the same 
as for previous observations (e.g. Bonfils et al. 2011b). Like 
GJ3634 -that hosts a super-Earth reported by our survey (Bonfils 
et al. 201 la)-, GJ 3470 is part of an extended sample of -300 M 
dwarfs specifically targeted to search for short-period planets, 
and the subset that transit. 

We collected 61 radial- velocity observations of GJ3470 be- 
tween 2008 Dec 08 and 2012 Jan 14. All are 900-s exposures 
except the first one that is a 300-s exposure used to verify the 
star is suitable for a planet search (i.e. not too active, not too fast 
a rotator, nor a spectroscopic binary). Their median uncertainty 



2 



X. Bonfils et al.: A sub-Neptune transiting the nearby M dwarf GJ3470 



Table 2. Modeled and inferred parameters for the GJ 3470 system. 





Unit 


Prior 




Posterior 






Stellar parameters 






Stellar mass, M+ 


[M Q ] 






0.541 +0.067 


Stellar radius, R+ 


[Ro] 






0.503 +0.063 






Planetary parameters 






Orbital period, P 


[day] 


Jeffreys (P min = 3.33, P max = 


= 3.34) 


3.33714 +0.00017 


Systemic velocity, y 


[km/s] 


Uniform (y min = +26, y max = 


= +27) 


26.51691 +0.00053 


Radial-velocity semi-amplitude, K 


[m/s] 


Jeffreys (K Mn = 0A,K max = 


: 100) 


0.00901 +0.00075 


Orbital eccentricity, e 


[] 


Uniform (e min = 0, e max = 


0.1) 


< 0.051 (1-cr upper limit) 


Argument of periastron, oj 


[rad] 


Uniform O min = 0, o) max = 


--In) 


- 27r (unconstrained) 


Planet- to- star radius ratio, R p /R+ 


[] 


Uniform (min - 0, max - 


0.1) 


0.0755 +0.0031 


Transit depth, D 


[mmag] 






5.69 +0.47 


Scaled semi-major axis, a/R^ 


[] 


Uniform (min = 0, max = 


100) 


14.9 +1.2 


Orbital inclination, i 


[rad] 


Isotropic spin-orbit (i^n = 0, i max = tt/2) 


> 88.8° (l-o- lower limit) 


semi-major axis, a 


[AU] 






0.0348 +0.0014 


Transit times 










Mid transit, T tr 


[day] 






55953.6645+0.0034 


First contact, t\ 


[day] 






55956.9652+0.0029 


Second contact, f n 


[day] 






55956.9705+0.0028 


Third contact, t m 


[day] 






55957.0325+0.0051 


Fourth contact, t w 


[day] 






55957.0325+0.0051 


Planetary mass, M^ 


[M e ] 






14.0 +1.7 


Planetary radius, 


[Re] 






4.2 +0.6 


Impact parameter, b 


[] 






0.22 +0.16 


Planetary density, p p 


[g/cm 3 ] 






1.07 +0.43 


Planetary surface gravity, g p 


[m/s 2 ] 






7.9+1.8 


Equilibrium temperature 










T eq , A=0 


[K] 






615+16 


T^, A=0.75 


[K] 






435 +12 



Radial- velocity time series additive noise 
for Harps, 6 rv [m/s] 

Photometric time series additive noise 

for Trappist (2012, Feb. 26th), 6 ph ,i [mmag] 
for Trappist (2012, Mar. 7th), 6 ph ,2 [mmag] 
for Euler (2012, Feb. 26th), 6 ph , 3 [mmag] 
for Nites (2012, Apr. 12th), 6 pM [mmag] 

Baseline model (y = a + br + cr 2 , where r is, for 
for Trappist (2012, Feb. 26th), a 
,b 
, c 

for Trappist (2012, Mar. 7th), a 
,b 
, c 

for Euler (2012, Feb. 26th), a 
,b 
» c 

for Nites (2012, Apr. 12th), a 

, b 



Data parameters 



Modif. Jeffreys (£ min 


= 0,6 


— 0.1,6 max 


= 1000) 


< 0.40 (1-cr upper limit) 


Modif. Jeffreys {e^ n 


= 0,e 


= 0.1, 6 min 


= 1000) 


2.18 +0.19 


Modif. Jeffreys {e^ 


= 0,6 


= 0.1, 6 min 


= 1000) 


3.51 +0.18 


Modif. Jeffreys {e^ n 


= 0,e 


= 0.1, 6 min 


= 1000) 


1.33 +0.12 


Modif. Jeffreys {e^ n 


= 0,6 


= 0.1, 6 min 


= 1000) 


0.04 +0.04 



a given time series, the time elapsed since the 1st exposure) 

1.000100 ±0.000461 
-0.011+0.012 
0.097 + 0.068 

0.999810 + 0.000502 
-0.003 + 0.011 
0.049 + 0.053 

0.998510 + 0.000355 
0.068 + 0.010 
-0.487 + 0.061 

0.186031 +0.000083 
-0.0058 + 0.0027 
0.302 + 0.018 



is (Ti = 3.7 m/s and result predominantly from photon noise. 
We note this uncertainty agrees well with that of other M dwarfs 
with similar brightness. The GJ3470 time series (shown in the 
top panel of Fig. 1 and available electronically) has an observed 
dispersion cr e = 1.1 m/s, indicating intrinsic variability. The pe- 
riodicity of that signal is seen in a generalized Lomb-Scargle pe- 
riodogram (Fig. 1, bottom panel - Press et al. 1992; Zechmeister 
& Kiirster 2009) where prominent power excesses are exhibited 
at periods P ~ 3.33 and 1.42 day (alias of each other for a 1-day 
time sampling). The periodogram is shown in Fig. I, together 
with false-alarm probability levels computed with bootstrap ran- 
domization of the original data (see Bonfils et al. 2011b,a, for 
more detailed descriptions). 



To model the data with a Keplerian orbit we chose a 
Bayesian framework and employed a Markov Chain Monte 
Carlo algorithm (MCMC - e.g. Gregory 2005, 2007; Ford 2005). 
MCMC algorithms sample the joint probability distribution for 
the model parameters by evolving a solution (i.e. a set of param- 
eter values) in the manner of a random walk. At each step, a new 
solution is proposed to replace the previous solution. The new 
solution is accepted following a pseudo-random process that de- 
pends on the x 1 difference between both solutions such that so- 
lutions with a higher likelihood are accepted more often. Step- 
by-step, accepted solutions build a chain which, after enough 
iterations, reaches a stationary state. One can then discard the 
first iterations and keep only the stationary part of the chain. The 
distributions of parameter values of all the remaining chain links 
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then correspond to the targeted joint probability distribution for 
the model parameters. Our implementation follows closely the 
one of Gregory (2007) with several (10 in our case) chains run- 
ning in parallel. Each chain is attributed a parameter J3 that scales 
the likelihood such that chains with a lower ft value present a 
higher acceptance probability. We also paused the MCMC it- 
eration after each 10 steps and proposed the chains to permute 
their solutions (which was again accepted pseudo-randomly and 
according to the likelihood difference between solutions). This 
approach is reminiscent of simulated annealing algorithms and 
permit evasion outside of local minima and better exploration 
of the wide parameter space. Only the chain with J3 = 1 cor- 
responds to the targeted probability distribution. Eventually, we 
thus discarded all chains but the one with ft = 1. We adopted 
the median of the posterior distributions for the optimal parame- 
ter values, and the 68% centered interval for their uncertainties. 
Hence, the orbital parameters were notably P = 3.3368 ± 0.0004 
day, K = 9.2 ± 0.8 m/s, e = 0.16 ± 0.08 and inferior conjunction 
T tY = 2455983.42 ± 0.14 JD, modulo integer multiples of the 
period P. Solutions with period P ~ 1 .42 day were also found, 
albeit with a 10 times lower occurence. 

For a M+ = 0.45+0.1 1 M© star, the optimal parameter values 
correspond to a companion with a mass m p sin / = 12 + 2 M®. 



4. Photometric detection 

From the radial velocity orbit, we could predict hypothetic tran- 
sit times with a 0.14-day accuracy. From mass-radius relation- 
ships (e.g. Fortney et al. 2007), we also estimated the transit 
depth caused by the eclipse of a pure iron planet would be of 
the order of 1 milli-magnitude and measurable from the ground. 
We aimed to perform a photometric search for a +2-<x window, 
i.e. 95% of the posterior density function (PDF) of transit times. 

4.1. First transit ingress detection with Trappist 

We initiated the photometric search with Trappist and ob- 
served on 2012, February 26th between 00 h 25 m and 04 h 36 m 
UT. Trappist is a robotic 60-cm telescope installed at La Silla 
Observatory. It is equipped with a 2048 x 2048 15/mi CCD cam- 
era, provides a 22'x22' filed-of-view with a 0.65" pixel scale 
(Jehin et al. 2011). We used a z'-Gunn filter, made 764 expo- 
sures of 10 s each and recorded flux for GJ3470 together with 
carefully chosen comparison stars in the field. We applied the 
reduction procedure described in Gillon et al. (2011) to derive 
the differential photometry (which time series are available in 
ascii form electronically). We note that changing the compari- 
son stars did not affect the results of our reduction. In Fig. 2 
(top) we show the time series after applying binning by 0.005 
day (7.2 min) and after computing uncertainties by measuring 
the dispersion in each bin and divided by the square-root of the 
number of points. 

Right after the first night of follow-up we identified a drop 
in flux compatible with the ingress of a ~ 4-R® planet, seen in 
Fig. 2 to start at BJD=2,455,983.66 day. The transit seems to 
occur somewhat less than 2cr later than predicted. We inspected 
Digital Sky Survey (DSS) images to look for a possible back- 
ground star. The proper motion of GJ3470 is fast and for its 
present position, the 60-year old DSS images can confidently 
exclude stars more than 5 mag fainter than GJ3470. We noted 
the event is not compatible with the aliased 1 .4-d period and we 
pursued the follow-up with Trappist, the Euler-Swiss and Nites 
telescopes with adjusted ephemeris. 



4.2. Confirmation with EulerCam, Trappist and Nites 

From La Silla Observatory, and after our first (fortunate) night 
of follow-up, only one transit event could be visible before 
the end of the current season. The ingress was foreseen at 
BJD^2,455,993.68 (2012 March 7th, 04 h 19 m UT) at an airmass 
greater than 1.8. To secure the confirmation we made use of 
EulerCam in complement of Trappist. EulerCam is a 4kx4k 
CCD camera mounted on the Euler Swiss telescope that is also 
installed in La Silla Observatory. The field-of-view of EulerCam 
is somewhat smaller than the one of Trappist (15.7'x15.7') but 
the twice as large telescope aperture mitigates the scintillation 
by a factor ~1.5. 

Trappist started observing at BJD=2,455,993.509 
(00 h 13 m UT), recorded 943 10-s exposures, and stopped 
at BJD=2,455,993.743 (05 h 50 m UT) when GJ3470 reached 
an airmass z = 3.55. EulerCam started at BJD=2,455,993.559 
(01 h 25 m UT), recorded 223 50-s exposures, and stopped at 
BJD=2,455,993.733 (05 h 36 m UT) when GJ3470 reached an 
airmass z = 3.04. We used z'-Gunn filter with both telescopes. 
The reduction procedure used for the EulerCam is described 
in Lendl et al. (2012). The time- series photometry presented 
has been obtained with relative aperture photometry, where 
apertures and references were chosen carefully. As for the 
previous light-curve, both time-series were binned by 0.005 
day and shown in Fig. 2 (second and third curves, from top to 
bottom). 

On both light curves, we identified a 6.6+0.4 mmag drop 
in flux consistent with the transit of a ~4-/? e planet, with an 
ingress timing in sync with our predictions. Unfortunately, the 
star reached a too high airmass for photometric observations be- 
fore we could record the transit egress. 

To record a full transit event we continued the follow-up 
in the Northern hemisphere, at La Palma, where we predicted 
the visibility of few more events before the star would go be- 
hind the Sun. The La Palma observations were obtained with 
the 0.4m Nites telescope (McCormac 2012). Images were ob- 
tained in white light using a 1024x1024 deep depleted CCD 
giving 0.7"/pixel and a 12'xl2' field of view. The entire tran- 
sit was observed without filter and at high airmass. The data 
were reduced with standard IRAF routines. Observations started 
on BJD=2,460,030.357 (2012 April 12th, 20 h 34 m UT) and 
stopped on BJD=2,460,030.500 (2012 Apr 13th, 00 h 00 m UT). 
We recorded 586 exposures of 20 s which aperture photometry 
is also shown at the bottom of Fig. 2, after binning by 0.005 day. 
The light curve appears with a trend of large amplitude which 
is most probably caused by the color-difference between our tar- 
get and the comparison stars. The airmass-extinction depends 
on the stars' colors and our target is sensibly redder than the av- 
erage comparison star in the field, an effect that is difficult to 
mitigate for unfiltered observations. Nevertheless the light curve 
does confirm the transit event and, more importantly, does con- 
strain its duration. 



5. Joined RV+photometry modeling 

To measure physical and orbital parameters of the GJ 3470 plan- 
etary system, we pooled together the photometry and radial- 
velocity time series. As for the spectroscopic orbit (Sect. 3), we 
chose a Bayesian framework and used a MCMC algorithm. We 
modeled the data with a planet on a Keplerian motion around 
the star, with 1 8 variables for the parametrization : the systemic 
velocity y, the orbital period P, the RV semi- amplitude K, the 
eccentricity e, the argument of periastron qj, the time of passage 
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Fig. 2. Light curves of GJ 3470. Left panel corresponds to raw light curves (with arbitrary offsets for clarity). Right panel corresponds to de-trended 
light curves. From top to bottom, light curves were obtained with Trappist (on 2012 February 26th), Trappist (on 2012 March 7th), EulerCam (on 
2012 March 7th) and Nites (on 2012 April 12th). The optimal model is over-plotted with a black curve (see Sect. 5). 



at periastron Tq, the planet- to-star radius ratio R p IR+, the scaled 
semi-major axis a/R+ and the orbital inclination on the sky i. 
In addition, to allow for quadratic baselines in the photometric 
time- series we add 4x3 parameters, and to model the additive 
error (i.e. quadratically in excess of photon noise) we used 5 ad 
hoc variables, e rv , £ P h,i, £ p h,2, £ph,3 and 6 p h,4, for each one of the 5 
RV and photometric time series. 

As dictated by Bayesian statistics we gave a priori descrip- 
tions for the probability distributions of parameter values. We 
chose uniform probability distributions for y, e, oj, Tq, a/R+, 
R p /R+, log(P) and log(K) (so called Jeffreys priors for P and 
K). We chose i such that orbits have isotropic distributions. And, 
we chose distributions that are uniform in a logarithmic scale 
above a given threshold and uniform in a linear scale below (so 
called modified Jeffreys priors), for all 4 additive noise param- 
eters. With each prior distribution we defined a range of autho- 
rized values (labeled by min and max subscripts) and, a threshold 
value for the modified Jeffreys priors (labeled with subscripts). 
Only for the baselines parameters, that are all linear, we solved 



for the best x 1 analytically. We report our choices for the prior 
distributions, parameters ranges and threshold values in Table 2. 

In addition, we assumed a stellar mass-radius relationship 
by interpolating within the Baraffe et al. (1998) models grid 
(with [M/H]=0 dex, Y=0.275, Lmix=1.0 and Age=l Gyr), and 
considered the stellar mass determination from Sect. 2 as one 
more piece of observational data rather than a parameter to fit. 
Note that we chose not to include a description of the Rossiter- 
McLaughlin effect. It is thought to be small compared to the pho- 
ton noise and, more practically, no radial-velocity happened to 
be taken during a transit event. 

To enhance the convergence efficiency, we restricted the in- 
tervals explored for the orbital period. This allowed us to con- 
verge without using several chains of different temperatures run- 
ning in parallel. We checked that the interval we chose was much 
larger than the posterior distributions for orbital periods and, 
therefore, did not affect our results. The MCMC chain converged 
on a stationary solution after -10,000 steps. We continued iter- 
ations for 500,000 more steps and inferred the posterior distri- 
butions of parameters from those last steps only. We adopted the 
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Fig. 5. Diagram with measured masses and radii for known exoplanets. Blue empty circles are for Venus, Earth, Uranus and Neptune (from left to 
right). Purple points are for previously known transiting exoplanets detected from the ground and green points are for planets detected from space 
by Kepler or CoRoT. Small vertical ticks on the top of the figure help label each detection. Kepler, CoRoT and HAT planets have short labels, K, C 
and H, respectively. GJ 3470 b is over plotted with a red diamond. The gray curves show mass-radius relations for water, rock (MgSi0 3 perovskite) 
and iron, from top to bottom, respectively (Fortney et al. 2007). 



median of the posterior distributions for the optimal parameter 
values, and the 68% centered interval for their ("1 <x") uncertain- 
ties (Table 2). 

We also picture the optimal model emerging from our 
stochastic fitting with Fig. 3, where it is plotted over the phase- 
folded RV and photometry data. Fig. 4 moreover shows the pos- 
terior distributions for a subset of system parameters. The top 
panel of that figure shows that our estimate for the planet's im- 
pact parameter is correlated to the stellar radius. For a small star, 
only a central transit matches the duration of La Palma's event. 
A larger star allows for a wider range of impact parameters but 
too large a star is rendered improbable by our prior input on M* . 
Our global modeling also attributes a mass, radius and thus den- 
sity to the planet, M p = 14.0 ± 1.7 M®, R p = 4.2 ± 0.6R e and 
p p = 1.07 ± 0.43 g/cm 3 , respectively. 

We repeated the above analysis with polynomials of higher 
degrees for the photometric baselines. NITES photometric trend 
is the strongest and we thus paid a particular attention to the 
egress timings. We tested 2nd, 3rd and 4th order polynomials 
and found slight changes in the timings determinations (consis- 
tent with the 1 - cr uncertainties quoted in Table 2). We also 
found the stellar radius determination and correlated parame- 
ters M+, M p and R p are sensitive to the choice of baseline. 
With higher degree polynomials, smaller values are found for 
all parameters (e.g. M* = 0.46 ± 0.08M o , = 0.43 ± 0.07R o , 
M p = 12.0 ± 1.8M© and R p = 3.49 ± 0.72R e for the 4th-order 
polynomial). 



Both Fig. 4 and the susceptibility to baseline choices show 
well that most of the uncertainties for those planetary parameters 
are bounded to the large uncertainties on the stellar properties 
(radius or mass). Fortunately, the possibility to refine both stel- 
lar and planetary parameters with more precise light curves is a 
virtue of planetary transits (Seager & Mallen-Ornelas 2003). For 
instance, a light curve of higher quality, such as those produced 
with Spitzer, shall refine ingress and egress durations. This will 
result in a strong constraint on the impact parameter, resolve 
the b-R+ degeneracy and improve the precision of the planetary 
mass, radius and density. 

6. Conclusion 

We have presented the transiting planet GJ 3470 b discovered 
with Harps radial velocities and subsequent photometric follow- 
up with Trappist, the Euler-Swiss and the Nites telescopes. The 
planet detection adds to the small subset of low-mass planets 
(M p < 30 Me) with measured masses and radii (Fig. 5). Today, 
its bulk properties remain largely unconstrained, mostly because 
we have limited knowledge on the stellar properties. With the 
current mass and radii determinations, GJ 3470 b would be an 
ice giant and seems comparable to Uranus in our Solar System. 
Among transiting exoplanets, it is intermediate to, on the one 
hand GJ436b (Butler et al. 2004; Gillon et al. 2007b,a), HAT-P- 
26 b (Hartman et al. 2011) and HAT-P-11 b (Bakos et al. 2010) 
and, on the other hand, GJ 1214 (Charbonneau et al. 2009) and 
55 Cnc e (Winn et al. 2011; Demory et al. 2011; Gillon et al. 
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2012). It is actually close to the nominal mass and radius of 
Kepler-llc (Lissauer et al. 2011b) and Kepler-20c (Gautier 
et al. 2012) which have loose mass determinations. 

Resolving the parameter degeneracy is likely to place 
GJ 3470 b as a remarkable planet. Indeed, on the one hand, if 
future observations attribute a small radius (e.g. ~ 0.45 R ) 
to GJ3470, the planet will be of low-mass (M p ~ 12 M©) and 
small (R p ~ 3.8 Re), in a yet unpopulated mass-radius domain. 
On the other hand, if future observations confirm or inflate the 
moderate radius of GJ 3470 (7?* > 0.50 R©), the planet will have 
an unusual low density (p p < 1.07 g/cm 3 ) compared to other 
ice giants (e.g. p = 1.3 and 1.6 g/cm 3 for Uranus and Neptune, 
respectively). 

Perhaps more importantly, the GJ3470 system has fa- 
vorable attributes for follow-up characterization (V mag =12.3, 
K mag =7.99, D=0.62%). Among super-Earths and exo-Neptunes 
only 55 Cnce, GJ436b and HAT-P-11 b transit stars brighter in 
both V- and K-bands. Among those 3, only GJ 43 6 b produces a 
larger transit depth (see Fig. 6). 

We thus encourage follow-up measurements of GJ 3470 b. 
Improved transit light curves should increase the precision of 
the planetary and stellar parameters. Chromatic light-curves (e.g. 
Bean et al. 2011; Berta et al. 2012) could probe the physico- 
chemistry of GJ3470b's atmosphere. On longer time-scales 
light curves could also measure GJ 3470's rotation and constrain 
its age with gyrochronology relations (Barnes 2007; Delorme 
et al. 2011). Astrometry could give a trigonometric parallax and 
improve knowledge of GJ 3470's distance. Finally, as the orbits 
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Fig. 4. Posterior distributions for selected system parameters. Iso- 
contours delineate 68.2, 95.4 and 99.7% confidence intervals (red, or- 
ange and brown, respectively). The median, 68% and 99.7% confi- 
dence intervals of marginalized single parameters are also reported with 
continuous-black, dash-red and dot-yellow lines, respectively 



of multi-planet systems are thought to have low mutual inclina- 
tions (Lissauer et al. 2011a; Figueira et al. 2012), the low value 
of planet b's impact parameter (< 0. 17 ± 0. 13) makes of GJ 3470 
an attractive target to search for additional transiting/cooler plan- 
ets. 
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